Experimental evidence for magnet or ot at ional instability in a helical magnetic field 
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A recent paper [R. Hollerbach and G. Riidiger, Phys. Rev. Lett. 95, 124501 (2005)] has shown 
that the threshold for the onset of the magnetorotational instability (MRI) in a Taylor- Couette 
flow is dramatically reduced if both axial and azimuthal magnetic fields are imposed. In agreement 
with this prediction, we present results of a Taylor- Couette experiment with the liquid metal alloy 
GalnSn, showing evidence for the existence of the MRI at Reynolds numbers of order 1000 and 
Hartmann numbers of order 10. 



The role of magnetic fields in the cosmos is two-fold: 
First, planetary, stellar and galactic fields are a product 
of the homogeneous dynamo effect in electrically con- 
ducting fluids. Second, magnetic fields are also believed 
to play an active role in cosmic structure formation, by 
enabling outward transport of angular momentum in ac- 
cretion disks via the magnetorotational instability (MRI) 
^ . Considerable theoretical and computational progress 
has been made in understanding both processes. The 
dynamo effect has even been verified experimentally, in 
large-scale liquid sodium facilities in Riga and Karlsruhe, 
and continues to be studied in laboratories around the 
world 2]. In contrast, obtaining the MRI experimentally 
has been less successful thus far Q. (p| claim to have 
observed it, but their background state was already fully 
turbulent, thereby defeating the original idea that the 
MRI would destabilize an otherwise stable flow.) 

If only an axial magnetic field is externally applied, 
the azimuthal field that is necessary for the occurrence 
of the MRI must be produced by induction effects, which 
are proportional to the magnetic Reynolds number (Rm) 
of the flow. But why not substitute this induction process 
simply by externally applying an azimuthal magnetic field 
as well? This question was at the heart of the paper 
where it was shown that the MRI is then possible with 
far smaller Reynolds (Re) and Hartmann (Ha) numbers. 
In this paper we report experimental verification of this 
idea, presenting evidence of the MRI in a liquid metal 
Taylor- Couette (TC) flow with externally imposed axial 
and azimuthal (i.e., helical) magnetic fields. 

The heart of our facility 'PROMISE' (Potsdam 
i?Ossendorf Magnetic /notability Experiment) is a 
cylindrical vessel made of copper (see Fig. 1). The use 
of copper was motivated by the fact that the instability 
usually occurs at lower Reynolds and Hartmann numbers 
for the case of ideally conducting boundaries than for 
non-conducting boundaries Q- The inner wall is 10 mm 
thick, and extends in radius from 22 to 32 mm; the outer 
wall is 15 mm thick, extending from 80 to 95 mm in r. 




FIG. 1: Sketch (left) and photograph (right) of the central 
module of the PROMISE facility. V - Copper vessel, I - Inner 
cylinder, G - GalnSn, U - Two ultrasonic transducers, P - 
Plexiglass lid, T - High precision turntables, M - Motors, F 
- Frame, C - Coil, R - Central rod. The dimensions are in 
mm. When the experiment is running, the rod R in the right 
picture goes through the center of the cylinders. 



This vessel is filled with the alloy Ga 67 In 20 5 Sn 12 - 5 , which 
is liquid at room temperatures. The physical properties 
of GalnSn at 25 °C are: density p = 6.36 x 10 3 kg/m 3 , 
kinematic viscosity v = 3.40 x 10 -7 m 2 /s, electrical con- 
ductivity a = 3.27 x 10 6 (Q m)" 1 . The magnetic Prandtl 
number is then Prn = \±§av = 1.40 x 10 -6 . 
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FIG. 2: Measured axial velocity in dependence on time and 
depth for different coil currents. 



The copper vessel is fixed, via an aluminum spacer D, 



on a precision turntable T; the outer wall of the vessel 
thus serves as the outer cylinder of the TC flow. The in- 
ner cylinder of the TC flow is fixed to an upper turntable, 
and is immersed into the GalnSn from above. It has a 
thickness of 4 mm, extending in radius from 36 to 40 
mm. There is thus a gap of 4 mm between this immersed 
cylinder and the inner wall of the containment vessel. 
The actual TC flow then extends between ri n = 40 mm 
and r out = 80 mm. At the moment, the upper endplate 
is a plexiglass lid P fixed to the frame F. In contrast, the 
bottom is simply part of the copper vessel, and hence 
rotates with the outer cylinder. There is thus a clear 
asymmetry in the endplates, with respect to both their 
rotation rates and electrical conductivities. 

Note also that the accuracy of this setup is not quite 
at the ~ 10 -2 mm level that can be achieved in ordi- 
nary, hydrodynamic TC experiments, e.g. 7]. For ex- 
ample, in order to ensure a well-defined electrical con- 
tact between the GalnSn and the walls, it is necessary 
to intensively rub the fluid into the copper. The preci- 
sion is then more like 10 _1 mm. The Reynolds number 
Re = 27r/i n rin(r ut — f\ n )/v is also O(10 3 ), considerably 
greater than the Re c = 68 value obtained in the nonmag- 
netic problem (with stationary outer cylinder). 

Axial magnetic fields of order 10 mT are produced by 
a double-layer coil (C). The omission of windings at two 
symmetric positions close to the middle resulted from 
a coil optimization to maximize the homogeneity of the 
axial field throughout the volume occupied by the liquid. 
This coil is fed by a power supply that can deliver up 
to 200 A. The azimuthal field, also of order 10 mT, is 
generated by a current through a water-cooled copper 
rod R of radius 15 mm. The power supply for this axial 
current delivers up to 8000 A. 

At present, the measuring instrumentation consists ex- 
clusively of two ultrasonic transducers with a working fre- 
quency of 4 MHz; these are fixed into the plexiglass lid, 
15 mm away from the outer copper walls, flash mounted 
at the interface to the GalnSn. They provide full axial 
velocity profiles along the ultrasound beamline by means 
of Ultrasound Doppler Velocimetry (UDV) |S| . A similar 
measuring system for the axial velocity components in 
TC experiments was first described by Takeda 9]. The 
comparison of the two signals from two opposite sensors is 
important in order to clearly distinguish between the ex- 
pected axisymmetric (m = 0) MRI mode Q and certain 
m = 1 modes which also play a role in some parameter 
regions of the experiment. 

Typically, the duration of the experimental runs was 
1900 sec, after a waiting time of one hour. Such a long 
waiting time was chosen not only due to the hydrody- 
namic gap diffusion time, but also due to our numerical 
predictions of rather small growth rates of the MRI mode 
in helical fields (see also jlflj). 

One of the features of TC flow under the influence of 
helical magnetic fields is the replacement of the Taylor 
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FIG. 3: Power spectra of the v z fluctuations, averaged over 
the depths between 220 and 308 mm, for 5 different coil cur- 
rents. The frequencies of the inner and outer cylinder are 
indicated. The MRI is visible at 50 and 75 A. In addition to 
that, an m — 1 mode is also present. 



vortex flow by a travelling wave 0, Q . This travelling 
wave appears already at fi := /out //in = 0, although with 
a very low frequency. With increasing /i, the frequency 
increases and typically reaches a value of 0.2/i n at the 
Rayleigh line /i Ray = (r in /r out ) 2 =: rj 2 (we have rj = 0.5 
here, hence /iR ay = 0.25). The crucial point now is that 
under the influence of helical magnetic fields, the critical 
Reynolds numbers remain relatively small increasingly 
far beyond the Rayleigh line, where nothing special hap- 
pens. Typically, this shift to larger \i becomes larger for 
increasing values of the ratio (3 := B^{r = r- in )/B z of 
azimuthal field to axial field. 

Another typical feature of the MRI is that, for fixed 
Re, it sets in at a certain critical value of the Hartmann 
number Ha = B z (ri n (r out — r^cr/pv) 1 ^ 2 , and disappears 
again at some higher value. In this letter we focus ex- 
clusively on experimental results which substantiate this 
behavior. 

All results presented in the following are for rotation 
rates of f- in = 0.06 s _1 and / ou t = 0.0162 s _1 , i.e. for 
H = 0.27, beyond the Rayleigh line /iR ay = 0.25. Figure 
2 documents a selection of seven experimental runs for 
coil currents J co ii between and 120 A. In each case, the 
axial current I vo & was fixed to 6000 A. The color coding 



of the plots indicates the (negative) axial velocity com- 
ponent measured along the ultrasound beam, from which 
we have subtracted the depth-dependent time average in 
order to filter out the two Ekman vortices which appear 
already without any magnetic field. These vortices, char- 
acterized by inward radial flows close to the upper and 
lower endplates, show up in our UDV data in the form 
of positive and negative axial velocities, separated ap- 
proximately at mid-height by a rather sharp boundary. 
This boundary most likely corresponds to a jet-like radial 
outflow in the center of the cylinder, as discussed in [ll|. 

For J CO ii = 30 A we already observe a traveling wave- 
like structure which is, however, still restricted to the 
middle part of the cylinder. One might speculate that, 
due to the (jet-like) radial outward flow there, fluid with 
lower angular momentum is transported outward, which 
leads to a steeper decrease of and hence to a (locally) 
lower value \i than what would be expected from the 
rotation ratio of outer to inner cylinder |ll[. 

For increasing i C oii? this traveling wave becomes more 
and more dominant, until at 75 A it fills essentially the 
entire cylinder. Increasing 7 CO ii even further though, at 
100 A this wave ceases to exist, and we again have a 
rather featureless flow. This is shown more quantita- 
tively in Fig. 3, where we show the average over the 
depths 220-308 mm of the power spectral density for the 
axial velocity for 5 selected values of 7 CO ii. A feature com- 
mon to all five curves is the appearance of the rotation 
rates of the inner and outer cylinders. This reflects cer- 
tain geometric imperfections of the facility, probably in 
the form of cylinder eccentricities or metal oxides sticking 
to some parts of the walls. There is also a certain peak 
approximately at the mean of the inner and outer cylin- 
der frequencies, which, on closer inspection of the two 
transducer signals, turns out to be a non-axisymmetric 
m = 1 mode. Most interesting for us here, however, is 
the MRI mode at f / f\ n ~ 0.1 — 0.2. These frequencies 
have been analyzed in detail for all values of 7 CO ii chosen 
in the experiment. 

We now interpret these experimental data in the con- 
text of numerical predictions, obtained and cross-checked 
by various independent codes 0, 0, ^| for the solu- 
tion of the linear eigenvalue problem in unbounded cylin- 
ders. First, Fig. 4a shows the dependence of the critical 
Reynolds number on the Hartmann number. We observe 
the typical MRI behavior, with the minimum value of 
Re cr [t occurring at Ha = 7.7. Fig. 4b compares the mea- 
sured frequencies of the travelling wave, normalized to 
the rotation rate of the inner cylinder, with the com- 
puted ones. The latter are actually given in two versions, 
at i?e cr it, and at the experimental value Re = 1775. It is 
interesting that the difference between the two curves is 
relatively small. Hence it might be anticipated that the 
normalized frequency in the nonlinear (saturated) regime 
is probably also not far from these values. The mea- 
sured frequencies are in reasonable correspondence with 
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FIG. 4: (a) Computed critical Reynolds numbers for varying 
/coii (and hence Hartmann numbers), for fixed 7 ro d = 6000 A. 
(b) Frequency /wave of the travelling MRI wave, normalized 
to the rotation rate fi n of the inner cylinder. The ratio of 
azimuthal field to axial field /3 is also shown. 



the computed ones, and show a similar behavior, with a 
maximum close to Ha =10. 

In summary, we have provided experimental evidence 
for the existence of the MRI in current-free helical mag- 
netic fields, by showing its appearance in a certain inter- 
val of Hartmann numbers, in reasonable agreement with 
numerical predictions. Further experimental results for 
other combinations of the governing parameters will be 
published elsewhere. Certainly, much work remains to 
be done, including a detailed investigation into the role 



of the magnetic axial boundary conditions. For later ex- 
periments, a symmetrization of these boundaries is en- 
visioned. Connected with this, a suppression of the Ek- 
man vortices by means of split rings (proposed in 
see also [l4|) may also help to avoid special effects in the 
mid-height of the cylinder. 
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